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Replacing the perturbative vacuum by the nonperturbative vacuum, we 
extend a recent development of a quantum field theoretic framework for scalar 
and pseudoscalar meson mixing. The unitary inequivalence of the Fock space 
of base (unmixed) eigenstates and the physical mixed eigenstates is inves- 
tigated and the flavor vacuum state structure is explicitly found. This is 
exploited to develop formulas for two flavor boson oscillations in systems of 
arbitrary boson occupation number. We apply these formulas to analyze the 
mixing of rj with rf and comment on the other meson-mixing systems. In 
addition, we consider the mixing of boson coherent states, which may have 
future applications in the construction of meson lasers. 
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I. INTRODUCTION 



The study of mixing transformations plays an important part in particle physics phe- 
nomenology |l| The Standard Model incorporates the mixing of fermion fields through the 
Kobayashi-Maskawa [0 mixing of 3 quark flavors, a generalization of the original Cabibbo 
H mixing matrix between the d and s quarks. In addition, neutrino mixing and oscillations 
are the likely resolution of the famous solar neutrino puzzle M\. In the boson sector, the 
mixing of K° with K° via weak currents provided the first evidence of CP violation || 
and the B°B° mixing plays an important role in determining the precise profile of a CKM 
unitary triangle |J in Wolfenstein parameter space fTj. The rj rf mixing in the SU(3) 
flavor group also provides a unique opportunity for testing QCD and the constituent quark 
model. Furthermore, the particle mixing relations for both the fermion and boson case are 
believed to be related to the condensate structure of the vacuum. The non-trivial nature of 
the vacuum is expected to hold the answer to many of the most salient questions regarding 
confinement and the symmetry breaking mechanism. 

The importance of the fermion mixing transformations has recently prompted a fun- 
damental examination of them from a quantum field theoretic perspective ||. A similar 
analysis in the bosonic sector has also been undertaken |J. However, more recent analysis 
|T0|j on the fermion mixing indicated that the previous result || based on the perturbative 
vacuum is only the approximation with respect to the exact one based on the nonpertur- 
bative (flavor) vacuum. In this work, we show that the same is true for the bosonic sector. 
Upon the completion of our work, we notice that the same conclusion was also drawn in a 
recent literature In our work, however, the orthogonality between mass and flavor vacua 
is shown in a straightforward algebraic method rather than solving a differential equation 



for the inner product of two vacua as presented in |TT] . As evidenced in the previous litera- 
tures P,|lOHl2ll the method of using a differential equation to prove the unitary inequivalence 
between the two Fock spaces has been known for some time and our algebraic method is a 
new development in this respect. Moreover, we analyze the structure of the nonperturba- 
tive flavor vacuum in a great detail contrasting to the fermion case. The details of flavor 
vacuum, its perturbative expansion in the mixing angle and also some clarifying remarks on 
the Green function method and the arbitrary mass parametrization are summarized in the 
accompanying Appendices. 

We begin in Section II with investigation of the vacuum structure using the relation 
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between the base eigenstate and the physical mixed-eigenstate fields. We derive the repre- 
sentation for Pontecorvo mixing transformation for boson case and explicitly calculate the 
flavor vacuum state structure in the quantum field theory (QFT). We then investigate the 
unitary inequivalence of the two Fock spaces - one is the space of mass eigenstates and the 
other is the space of flavor eigenstates. In Section III, the ladder operators are constructed 
in the mixed basis. These are used to derive time dependent oscillation formulas for 1-boson 
states, n-boson states, and boson coherent states. Consequences from the replacement of 
the perturbative vacuum by the exact nonperturbative (flavor) vacuum are demonstrated. 
Section IV is devoted to study specific cases in our formalism, such as the r] 7/ system. We 
show the numerical differences between the two results: one from the perturbative vacuum 
and the other from the nonperturbative vacuum. Conclusions and discussions follow in 
Section V. In Appendix A, we present a derivation of an explicit expression for the flavor 
vacuum operating the ladder operators of particle and antiparticle to the vacuum of mass 
eigenstates. In Appendix B, we discuss the region of validity for a perturbative expansion of 
the flavor vacuum. In Appendix C, we make some clarifying remarks on the Green function 
method and the arbitrary mass parametrization discussed in recent literatures . 



II. THE MIXING RELATION AND VACUUM STRUCTURE 



We start our analysis by considering the Pontecorvo mixing relationship [14[] for two 
fields: 

4> a = cos6ipi + sin 6^2 

4>/3 = — sin 6^1 + cos9ip2, (2.1) 

where y9i j2 are the free fields with definite masses mi ;2 and Qj/ g are the interacting fields 
with definite flavors a, ft, respectively. The above mentioned relationship naturally arises 
by considering the mixing problem for the two quantum fields with the lagrangian of the 
form 

L = L O)a + L OtP -X(^ + ^ a ), (2.2) 

where A), a (j9) are the free flavor- field lagrangians (i.e. Lo,a{p) — K^lrm^aGS) — 
m a(/3)^Q(/3)^a(/3))) an d ^ * s ^ ne coupling constant responsible for mixing. It is straightforward 
to show that the above lagrangian can be immediately diagonalized by the transformation 
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given by Eq. (|2.1|) with an appropriate choice of mixing angle 9 [10 . The parameters of 
diagonalized lagrangian can then be expressed in terms of flavor-field masses (m a ,m/3) and 
interaction constant (A); i.e. 

4A 

tan(20) = — 2 (2.3) 



and 



(m 2 a + m 2 j )±J (ml - m}) 2 + 16A 2 
™\ 2 = N — ■ (2-4) 

The free mass-fields tpx >2 can be written explicitly as usual: 

I2ti(k) 

where ap and bp are respectively the particle and anti-particle ladder operators for the free 
mass-fields and they satisfy the standard equal-time commutation relationships: 

hA^ = k^'- ( 2 - 6 ) 



Here, kx = k x — k ■ x and e»(/c) = k (k) = yk 2 + mf. For the spin-0 case, free mass-field 
amplitudes up and vp are just numbers, i.e. = vp = 1- The interacting flavor-fields 
given by Eq.( |2.1| ) are the solutions of the Euler-Lagrange equation for Eq.( |2.2| ) and therefore 
can be completely determined in terms of the two free spin-0 fields given by Eq.( |2.5| ) and 
the mixing angle 9. 

This, however, gives rise to highly nontrivial relationship between the Fock-space of 
free-fields and that of interacting fields. To build the Fock space of flavor-eigenstates we 
consider the representation of the transformation consistent with Eq. (|2.1| ) in the Fock space 
of mass-eigenstates. Using the Baker-Hausdorff lemma, we can write the generator of this 
transformation as 

S = J d?x(4>i(x)<pp(x) + (j) ] p(x)$ a (x) - 4>l(x)cj) a (x) - (j) j a (x)<pp(x)) 

= f d 3 x((p\(x)ip 2 {x) + ipl(x)ipi(x) - ipl(x)ipi(x) - ip\(x)ip 2 (x)) (2.7) 

so that Eq.( p.l[ ) can be written in the form: 
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<Pp = e~ iS9 ip 2 e l 



(2.8) 



with the transformation operator given by 

G{6) = e^ 9 . (2.9) 

The similarity transformation given by Eq. Q2.8| ) relates the free field operators <^ 12 to the 
interacting fields <p at p. These relationships can be obtained from the requirement of the 
inner product conservation 

< a \(p\(3 >=< a\ip\b >, (2.10) 

where a linear transformation of state vector according to 

\a >= G~ x {9)\a > (2.11) 

relates the two Hilbert spaces: i.e. mass-eigenstate space TLip and flavor-eigenstate space 
7~la,/3 — G~ 1 (9)T-Ci t 2- The transformation given by Eq.( |2.8| ) can also be viewed as the "rota- 
tion" of basis in the Hilbert space of quantum states diagonalizing the bilinear lagrangian 
given by Eq. fl2.2| ). 

The operator S can then be written in terms of ladder operators a& and as follows: 

$ = E 5 {7+(a&4 2 + b -ki b lk2 ~ a li a k2 ~ b Ui b -k2) + 
% 

+ 7-(o&6-a + H2 b -ki ~ 4x b -k2 - b -ki a l)}> ( 2 - 12 ) 



where we denote 7+ = + ^ and 7 _ = ^ jgj - Here 7+ - 7- = 4. From 

Eq.( |2.12| ) we note that Eq. ( |2.8| ) makes each cluster fi^, defined by linear superposition of 
operators (a% v a^ 2 , b_~ , b ? ), transform into itself, i.e. e~ lSe Vt^e lSe = Fl^. The same can be 



given by Eq. fl2.8|) within each cluster with a specific momentum k: 



said about its hermitian conjugate Ql. This means that we can consider the transformation 



S% = I {7+(«fel4 2 + b -kl b lk2 ~ 4i a k2 ~ b Ul b -k2) + 

+ l-iHx b -k2 + H2 b -ki ~ 4i & U ~ fe -«4 2 )}- ( 2 - 13 ) 

Thus, the total transformation is given by e lSe = T7fc e lS i e . It is also convenient to express 
as := \/2(Tt + T^) with the operator defined by 

Tu = -^=( 7+ (4i^ - b -kA %2 ) + 7-(4i 6 -& - °& 6 -Si)). ( 2 - 14 ) 



where the commutation relation [T^,Tl] = 1 is satisfied between T k and TjJ just the same 
way as the particle creation and annihilation operators satisfy the commutation relations. 
With the operators T fc and T^, we can directly apply Eq.( |2.9|) to the mass-eigenstate vacuum 
and obtain 

|0; 9, k > a ,,= £ { ~ l6 f T ± C l n (9)(flyff l \0 > 1)2) (2.15) 



where C l n (6) are the generalized binomial coefficients that can be found after appropriate 
order ings of T and T' are carried out. In the expression given by Eq. fl2.15l) , one can treat 
the operators (T-!) as the annihilation (creation) operator of the vacuum fluctuation. 

For simplicity we now suppress the momentum notation in the ladder operators and take 
the flavor vacuum state in the most general form 

|0;fl>= Y, C nlmk (6)(a{) n (al)\bl 1 ) m (bU) k \0> . (2.16) 

n,Z,m,fc 

Applying the flavor annihilation operators to this vacuum, we obtain an infinite number of 
coupled linear equations for the coefficients C n i mk {9) and solve these equations in Appendix 
A. As shown in Appendix A, we find 

|0, 9 >= zY-^iZnalbl, + Z 12 a{bl 2 ) n (-Z n albU + Z t2 aW-i) l \0 >, (2.17) 

where Z =< 0\0,6 > is the normalization factor to be fixed by < 8;0\0;6 >= 1 and the 
coefficients Z n and Z 12 are given by 

7 + 7_ sin 2 9 7+7- sin 2 9 



Z 



li 



4(cos 2 9 + \ sin 2 9) 4(1 + \ sin 2 i 
— 7_ sin 29 ~1- si n 2$ 



Zj\2 — 2 — 2 • (2.18) 

4(cos 2 9 + \ sin 2 9) 4(1 + \ sin 2 9) 
Here, we note that the coefficients Zu and Z 12 (see Eq.( p.l8|) ) can be written as 

z u = x-x } z l2 = x-y, (2-19) 

where 

7_ sin 9 

X = 



2\1 



x 



4 

7+ sin 9 



7^ sin 2 ( 



2\jcos 2 9 + 

) 

(2.20) 



i 

cos 9 



cos 2 9 + 



7^ sin 2 ( 



Thus, the flavor vacuum state given by Eq. fl2.17|) can be rewritten as 



Ti-\-l 

\0 i e>=Z^2^ (zat&ti + ya{bl 2 ) n (-xa\bU + ya\bl^ |0 > . (2.21) 

71, 1 

This result can be further reduced as 

IM >= 2£^(aMr(«M)'|0 >, (2.22) 



, n!Z! 

n,l 



by defining new ladder operators: 



ci = + yb^ 2 , 

C2 = -xb-2 + yb-x, (2.23) 

where it is easy to check that [c\, c\) = x 2 + y 2 = 1, [c2, <4] = 1> [ c i; c 2] = 0. Now, it is 
possible to compute directly the value of Z, because 

y 2(n+Z) / \ 2 Z 2 

<0,9|0.«>=Z 2 i:^ ]r n!^=^(5:x 2 j =(Tf^ (2-24) 

and the flavor vacuum is normalized to be one,vi.e. < 0, 9\0, 9 >= 1. Thus, we find 

Z = l- X 2 = : - 2 l . 2fUA . (2.25) 
1 + 7_ sin 0/4 

We see that the flavor vacuum state indeed exists in the Fock space of mass-eigenstates and 
the normalization factor Z is finite but less than one for any value of 7_ sin 6 in the exact 
vacuum treatment. The same has been obtained in Ref. [IT] solving a differential equation 



of < 1 0, 9 >=< 0\G(6)\0 >. While such method of derivation using a differential equation 



has been known for some time pilOHT2|, our algebraic method presented in this work is a 
new development. 

This proves then the unitary inequivalence between the two Fock spaces of mass and 



flavor in the infinite volume limit following the procedure discussed in Ref. |TT 

V r , 

lim < O|O,0 > af s= lim exp( — - / d 3 k\nZ) = 0, (2.26) 
for any time. While we agree with Ref. [O] on the point that only the infinite volume limit 



can warrant the unitary inequivalence even in the boson case, we note that the perturbative 
expansion of the exact vacuum in the boson case is dramatically different from the case of 
fermion. Although the normalization factor Z is a finite function for all values of 7_ sin 9, 
we observe that this expression given by Eq. (|2.25 ) has singularity on the complex plane 



at 7_sin# = 2i. This is in a remarkable difference from the fermion case where the cor- 
responding result Zf ermion = 1 — 7^sin 2 #/4 doesn't have any singularity on the complex 
plane. Thus, the flavor vacuum |0, 9 > in terms of series in 9 shall have a critical point and 
this would result in the divergence of the Taylor series expansion for < 0|0,# > in powers 
of 9 because such expansion only makes sense for small 9 values. As we explicitly show in 
Appendix B, the series given by Eq.( |2.15D is indeed divergent in the region 7_# > 2. Such 
divergence doesn't occur in the fermion case. We also present some clarifying remarks in Ap- 
pendix C regarding on the Green function method and the aribitrary mass parametrization 



discussed in the previous literatures [10,11,13 



III. LADDER OPERATORS AND CONDENSATIONS 



In the previous section, we have built the representation of the mixing transformation 
given by Eq.( |2.1| ) in the operator space of ^1,2, where the action of mixing is given by the 
similarity transformation given by Eq.( [2.8|) . We also considered the representation defined 
by operating G~ 1 (9) in the Fock space and showed the unitary inequivalence between the 
two (mixed and unmixed) Fock spaces in the infinite volume limit. 

Let us now further investigate these representations to come up with physically measur- 
able quantities. The fields <pi$ are defined by a superposition of ladder operators a\ y2 and 
61,2 that form the basis in a linear Hilbert space of mass eigenstate fields. Using Eqs.( [2.5| ) 
and (|2.8|) , one can immediately obtain annihilation operators for the mixed (flavor) fields 
that are consistent with the Pontecorvo mixing relationship; 

a a ,p = G- 1 (6)a 1 , 2 G(e), 

b a ,p = G- 1 (9)b lt2 G(9). (3.1) 

This is also consistent with the definition of flavor vacuum as the lowest energy state, i.e. 

< 0, e\H(6)\0, 9 >=< 0\G(9)H(9)G- 1 (9)\0 >=< 0\H \0 >= 0, (3.2) 

where H(9) and Hq are the Hamiltonians of mixed-fields and unmixed-fields, respectively. 
Straightforward application of Baker-Hausdorff lemma to Eq. (|3.1|) yields: 

, sin9 t 
a a = a x cos 9 + — — (7 + a 2 + J-bL 2 ), 

sin ^/ ,t \ 

a/3 = a 2 cost* + — — (-7+ai + 7- & -i), 
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sin6> , t\ 

0-a = 0-1 COS 6" + —^—(7+0-2 + 7- a 2) ; 

b-p = 0- 2 cos0 + ^^(-7 + 6_! + 7-4). (3.3) 



It is also not difficult to reverse Eq.( |3.3| ) in order to obtain how the mass-eigenstate ladder 
operators are expressed in terms of flavor ones. Using the above relationships we can also 
find the time dependence of the flavor-eigenstate ladder operators in the Heisenberg picture 
since the time evolution of mass-eigenstate ladder operators are given by: 

a h2 (t) = e iflot a h2 e' lflot = e^'a^, 

b h2 (t) = e i&0 % 2 e- iAot = e~ iei ' 2t b 1>2 . (3.4) 

In particular, after introducing more compact notation 

sin # 7+ sin fl 7 _ 
C = cos 0; S + = ; d_ = , (3.5) 



we find 



a at = (C 2 e~ ieit + S 2 + e~ ie2t - S 2 _e it2t )a a + CS + (e~ ie2t - e~ ieit )ap+ 
S + S^{e ie2t - e~ i€2t )bl a + CS_(e it2t - e'^bl^ 

a pt = (C 2 e~ ie2t + S 2 + e~ ieit - S 2 _e Ult )ap + CS + (e~ le2t - e~ iei *)a a + 
S + S_(e- i£lt - e i£lt )bl p + CS^{e ieit - e~ it2t )bl a , 
b_ at = (C 2 e~ ieit + S 2 + e~ ie2t - Sle it2t )b_ a + CS + {e~ ie2t - e~ ieit )b^p+ 
S + S_(e ie2t - e~ i€2t )ai + CS_{e ie2t - e"" 1 *)^ 



-fit 



(C 2 e~ ie2t + S 2 + e~ ieit - S 2 _e ieit )b_ p + CS + ( e - lt2t - e- leit )b^ a + 

S + S_ (e- i€lt - e i€lt )al + CS_ {e i£lt - e _<eat )at , (3.6) 



from which we can also obtain the unequal-time commutation relationships: 

[a at a at] = [°— aj °— a t] — C e 1 + S^e 2 — S_e 2 = A aa , 
[a p , a pt ) = [b- P , bl pt ] = C 2 e^ 2t + S 2 + e^ - Ste'^ = App, 
[ap, al t } = [a a ,ap t ] = [b-p,b_ at ] 

= [b- a , blp] = CS + {e u * - e ieit ) = A Pa , 
\b-(3,a ai \ = [a[3,b- a t] = —[b- a ,apt]* 

= -[a a ,b^ t ]* = CS-(e Uat - e"* 1 ') = Ap a , 
[b. a ,a at ] = [a Q ,b. at ] = S + S.(e le2t - e~ te2t ) = A aa , 

[b-p, ap t ] = [ ap , b-pt) = S+S-ie-** - e^) = Ap p . (3.7) 
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All other commutators are either zeros or can be expressed in terms of the above ones. 
Eqs . (|3 . 3| ) , ( P . 6| ) , ( 377 ) in fact define all the dynamics of Pontecorvo mixing for two quantum 



fields. To show how these relationships can be used to calculate the dynamical parameters 
of the mixed (interacting) fields, one can consider the time evolution of cluster defined 
in Section II. As discussed in Section II, however, this cluster is invariant under the G~ 1 {9) 
transformation. Thus, we can consider Qj: with a particular k independently from all other 
momentum values. 

We now calculate the number of particles with a definite mass condensed in the flavor 
vacuum state |0' >= |0 > a ,p. Let's consider the condensation of the particle with a definite 
mass, for example Z\ =< 0'|iVi|0' >. Using the inverse relation of Eq. (|3.3| ): 



sin 9 + 

a ± = a a cos 9 — (7+6^ + 7_fr_ /3 ), (3.8) 

we can get 

• 2/1 2 

Z x =< 0'|4ai|0' >= Sm - 7 ~ < tflb-pblplQf > . (3.9) 

One can show that the same result is true for Z 2 =< 0'|iV 2 |0' >• Thus, the condensate 
density of particles with a definite mass in the flavor vacuum is given by: 

Zl = Z 2 = Si = (3.10) 

Apparently the condensate densities for particles with definite flavor in the mass vacuum, 
i.e. < 01^(^)10 >, are also given by S 2 ^. Let us now consider the number of particles with 
a definite flavor in the flavor vacuum, for example Z a (t) =< 0'\N a (t)\0' >. Using Eq. (|3.6| ), 
one can easily show that 

Z a (t) =< 0'\(S + S^(e te2t - e~ te2t )bl a + CS_(e te2t - e~ ieit )bl^j ■ 

(S + S4e ie2t - e~ lt2t )bl a + CS_{e l€2t - e^ 61 *)^)^' > (3.11) 

and thus 

Z a (t) = ASlS 2 + s\n\e 2 t) + ASlC 2 sin 2 (^^t). (3.12) 
Similarly, we get for the /3-particles 

Z p {t) = AS 2 _S 2 + sin 2 ( ei t) + ASlC 2 sin^^^t). (3.13) 
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We see that the number of particles with a definite flavor in the flavor vacuum is indeed not 
zero. This is due to the fact that flavor vacuum is not an energy eigenstate of the hamil- 
tonian H(8) and changes with the time translation producing and destroying coherently 
virtual particle/antiparticle pairs. It shows a significant difference from the ordinary quan- 
tum mechanical treatment without considerating the vacuum effect, which yields Z a tp\ = 
for any time. We emphasize that our flavor vacuum here is not perturbative but exact. This 
is different from the approach, where mass eigenstate vacuum |0 >i,2 is used instead of the 
flavor vacuum to generate a flavor eigenstate, e.g. \a >= aJjO >i,2- If the flavor vacuum 
|0' > was replaced by the mass vacuum |0 >i,2, then we would have obtained Z\ = Z 2 = 
instead of Eq. (|3.10|) . As discussed above in the exact vacuum treatment, the mass eigen- 
state vacuum is not annihilated by a a ^ operators. Indeed the term proportional to 0(7-) 
remains in the creation/annihilation operators, so that the accuracy in the order of 0(7^) 
can be expected from the results of exact vacuum approach compare to the perturbative 
vacuum approximation. The densities of vacuum condensation for antiparticles Z\^ Z 5 ^ 
are obtained same as the densities for the corresponding particles, i.e. Z lj2 , Z a ^ given by 
Eqs.(|3Tq),(FT2D and (|3T3| ). 

We now consider the flavor oscillations in time for a single particle with flavor a and 
momentum k. In the Heisenberg picture, the average number of particles with flavor a = a 
or (3 in the flavor state \ot >= aJjO' > is given by: 

< N a (t) >=< a\N a (t)\a >=< 0'\a a a{ t a at al\0' > . (3.14) 

In Eq. (|3.14j ), we note that we use the flavor vacuum to obtain exact result for the flavor 
oscillations. Later, in Section IV, we numerically compare the exact result with the previous 
approximate result 0. Using Eqs. ( [3 . 3|) , (|3 . 6| ) , (|3 . 7| ) , we directly apply the standard quantum 
field theoretic method. Since the flavor vacuum is annihilated by a Qj/ g, we move a a in 
Eq.( |3.14|) to the most right position and a) a to the most left position to annihilate the flavor 
vacuum. What is left is uniquely determined by the unequal time commutation relations 
given by Eq. ( |3.7| ) and we find 

< a\N at \a > = < 0'\aLa at \0' > +\\a a ,aL}\ 2 = Z a + \ A aa \ 2 \ 
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< a\N_ &t \a > = < / |6l at 6_ at |0' > +| [a a , b„ at ] | 2 = Z a + \A t 
< a\Npt\at > = < 0'\a' /3t ap t \0' > +|[a Q ,a^]| 2 = Zp+ \A Pa \ 2 ; 

< a\N_- m \a > = < / |6l jBt 6_ /9t |0 / > +|[a a ,6_^]| 2 = Z p +\Ap a \ 2 . (3.15) 

Using the notation of C, S±, our results are summarized as: 
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< a\N at \a > = 1 + 8C 2 S 2 _ sin 2 (^±^t) + 85^5j sin 2 (e 2 i) - 4C 2 S 2 sin 2 (^^ 2 -t), 



< a\N 0t \a > = AC 2 S 2 _ sin 2 (^±^t) + AS 2 S 2 + sm 2 (e 1 t) + 4C 2 S 2 sin 2 (^pt), 

< a\N_ at \a > = AC 2 S 2 _ sm 2 (^-^t) + 8S 2 _S 2 + sin 2 (e 2 t), 

< a\N_ 0t \a > = 8C 2 S 2 _ sm 2 (^^-t) + AS 2 _S 2 + sin 2 ( ei t). (3.16) 

As shown in Eg. ( |3.16| ),the time dependence of the average number of particles with a definite 
flavor is rather complicate. It contains oscillating contributions both from the a — > j3 
conversion and from the virtual pair creation in a dynamically "rotating" flavor vacuum. 
As discussed in Ref. M, the a — > (3 conversion process generates the term proportional to 
sm 2 ( ei ~ 62 t). The terms involving t\ + e 2 , ei, e 2 -frequencies in Eq.( p.l6| ) are, however, related 
to the creation of virtual pairs. For example, the virtual pair creation violates energy 
conservation within the uncertainty time, i.e. AEAt ~ 1 (in our units h — 1) and thus both 
creation and annihilation of, let's say, (a + p) virtual pair must occur within r pa ei ^ e2 time 
interval. Thus, the terms in Eq. ( p.!6|) involving e 1 + e 2 , ei, e 2 -frequencies can be related to 
the creation of different types of virtual pairs, while the terms involving e\ — e 2 is related to 
the actual a — > (3 conversion. 

Using Eq.( |3.16 ), we can also calculate the expectation value of the flavor charge operator 
defined by Q a ,/3 = N a)/3 - N_^,-p] 

<Q a > = 1 - AC 2 S 2 + sin 2 (^^) + ACS 2 _ sin 2 (^-±^), 

<Q p> = AC 2 S 2 + sm 2 (^^t)-AC 2 S 2 _s\n 2 (^-^t), (3.17) 

or with the conventional parameters, 

<Q a > =1-7+ sin 2 (20) sin 2 (^=-^) + 7 2 sin 2 (26) sin 2 (^-±^), 

<Qp> = 7 2 sin 2 (2fl) sin 2 (^^t) - 7 ! sin 2 (2#) sin 2 (^-^t). (3.18) 

From this result, one can also see that there is an additional term proportional to sin 2 ( 61 ^ <;2 t) 
to the usual Pontecorvo formula. As discussed above, the origin of this term can be under- 
stood as a contribution from the virtual pair creation in "rotating" vacuum. The correction 
term is of the order of 0( 7 2 ). As noted earlier, this may give a reason why it has been 
found neither in an ordinary quantum mechanical treatment nor in the approximate QFT 
treatment based on a perturbative vacuum. 

We also calculate the time evolution of coherent state for the two mixed quantum fields. 
Coherent state has the form 
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\Ca >= e Ca «\0' > . (3.19) 

Extending the above calculation for a single particle, it is not so difficult to verify that for 
the state containing n particles with flavor a can be given by 

< n\N at \n >=\< 0'\a n a N at (ai) n \0' >=< N at > +n\A aa \ 2 . (3.20) 

Besides n|A Qa | 2 which is simply n times the probability of a — > a transition, we see in 
Eq. ( |3.20|) that the condensate contribution is present adding the density of a particles from 
"rotating" vacuum. Applying this result directly to the coherent state expansion, we obtain 
the following expectation values of the number operator N^ a ^ t in the coherent state [Cat >: 



|2| a |2 
'aa | j 



< Ca\N at \Ca >= Z a + \C\ Z \A C 

< Ca\Np t \Ca >= Z p + |Cf |^V| 2 - (3.21) 

Thus, the expectation values of the flavor charge operator Q(a,p) — ^{a,p) — Nf_a,-i3) t urn 
out to be 



< Ca\Q a \Ca > = |Cf <Q a > 

= |C7| 2 (l-7?. 6 

< Ca\Qp\Ca > = \C\ 2 <Q P > 



Cf (1 - 7 ^ sin 2 (2^) S in 2 (^-^t) + 7 ! sin 2 (2^) s in 2 (^^t)), 



= |C| 2 ( 7 ^ sin 2 (2^) S in 2 (^^t) - 7 ! sin 2 (2^) s in 2 (^^t)). (3.22) 

As we can see in Eq. (|3.22| ), the vacuum contributions Zi a m are removed from the flavor 
charge expectation values and the results for the coherent state are simply \C\ 2 times the 
expectation values of the flavor charge for the single particle state. 



IV. APPLICATION TO REAL MESON STATES 



We now apply the results for time evolution of two mixing boson fields to the analysis of 
rj — rj' mixing system. The masses are taken to be 549 MeV and 958 MeV, respectively, and 
of course in the particle rest frame the energies in our formulas reduce to the masses. The 
phenomenologically allowed mixing angle (0su(S)) range of the r\rj system is given between 
— 10° and —23° fTSf , where the mixing angle 9su(3) is defined by Eq.(36) of Ref. fi~6f . This 
angle represents the breaking of the SU(3) symmetry, the eigenstates of which are already 
rotated —35.26° from uu + dd and ss to a = uu + dd — 2ss and (3 = uu + dd + ss. Thus, 
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our mixing angle is defined by 9 = 9 



SU(3) 



35.26°. Recent analysis of the rjr)' mixing angle 



using a constituent quark model based on the Fock states quantized on the light-front can be 
found in Ref. []TB and the references therein. The optimal value found for 9su(3) was around 
— 19° and thus 9 « —54°. We use these values in Eqs.( |3.16"D and ( |3.17| ) (or equivalently 



( to determine the evolution of definite flavor particle number and charge. 



"D 

O 
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FIG. 5.1. comparing population density evolution for k=0.1GeV 



30E-24 



In Fig. |5.1| , we present both < aliV^a > (thick solid line) and < Q a > (dotted line) as a 
function of time when the particle momentum is given by k = 0.1 GeV. For a comparison, 
we also show the previous approximate result (thin solid line) based on the perturbative 
vacuum || corresponding to these quantities noting that < a\N at \a > and < Q a > coincide 
each other in this approximation as one can see in Eqs.( |3.16j ) and ( |3.17[ ). As we shown in 
Fig JO, the population density < a\N at \a > (thick solid line) is completely distorted due 
to the interaction with the nonperturbative vacuum while the sinusoidal Pontecorvo result 
(thin solid line) is obtained for the approximate perturbative vacuum treatment. We see 
the large deviation up to 40% in < a\N at \a >. However, one cannot see the same level of 
deviation in < Q a > and the previous result || based on the perturbative vacuum seems to 
be a good approximation for the description of flavor charge oscillations modulo the accuracy 
of order 0(j^_). 
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FIG. 5.2. population density evolution for k=0.1GeV 



More details of our results on the time evolution of the particle number with the mo- 
mentum k = 0.1 GeV are shown in Fig. |5.2| , where the thick solid and dashed lines are 
< a\N at \a > and < a\Np t \a >, respectively, and the thin solid and dotted lines are respec- 
tively the antiparticle contributions of < a\Nat\a > and < a\Ng t \a >. The 7777' is one of 
the most severely mixed systems due to the great difference in masses of mixed particles. 
As we have stated earlier, the simple harmonic structure of average particle number usually 
obtained in quantum mechanics or in an approximate QFT treatment 15) is completely al- 
tered as a result of nontrivial interaction with the complicate vacuum. What we see is the 
superposition of two different cycles as described by Eq. (|3.16| ). From the initial moment of 
time the population of both a-particles (thick solid line) and /3-particles (thick dashed line) 
increases. Although the increase of number of /3-particles in system is well understood due 
to a — > (3 conversion, the initial increase of a-population is quite unexpected and caused 
by a — a production from vacuum. The contribution from this process however is rather 
fast so that the general tendency of exchanging between a and j3 particle states can also 
be seen quite well. In Fig. |5.2| , we also see the oscillations of the antiparticle number in the 
system. This effect is given in the order of 7^ and usually is absent in an approximate 
QFT treatment. This is entirely a QFT effect which cannot be obtained within the frame- 
work of quantum mechanics. In QFT, besides the beams of a and (3 particles moving in 
k direction, we necessarily have antiparticle beam traveling in the opposite direction. The 
population density in this beam is correlated with particle-beam so that the total flavor is 
preserved. The existence of beam is caused by " dynamically rotating" vacuum disturbance 



15 



at the initial time of a-particle emerging. One should also note that the existence of "recoil" 
antiparticle beam is preserved in the more general wave-packet QFT treatment of mixing 
problem. Thus, the mixed particle of definite flavor not only produces the usual oscillation 
of population density in time (or space |T7|]) but also is accompanied by emitting the beam 
of antiparticles traveling in the direction opposite to the beam of particles. These effects are 
in principle testable in the experiments. 
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FIG. 5.3. flavor charge oscillations 



In Fig. |5.3| , we also plot more details on the time dependence of flavor charge expectation 
value with the same momentum k = 0.1 GeV. The thick solid and dashed lines are < Q a > 
and < Qp >, respectively. One can see that they exhibit mainly the simple periodic structure 
similar to the approximate QFT-results and only slightly distorted due to interaction 
with vacuum. The amount of distortion is of 7^ order, i.e. about 10% for this case. 
Interesting feature is however presence of the regions where flavor charge of given sort of 
particles change sign which means that antiparticles outnumber the particles. The process 
can be physically understood as result of a — a production when number of a-particles is 
small due to a — > (3 transition. 
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FIG. 5.4. mixing amplitudes 

It is also interesting and experimentally testable that the efficiency of conversion processes 
and the flavor-vacuum disturbance depend essentially on the energy of original particle. The 
dependence is effective to the relativistic mass of particles so that the QFT-mixing effects 
are decreasing with the energy increase of flavor-particle. The distribution of intensity 
for simple quantum mechanical mixing and QFT-mixing is given by the relationship of 
amplitudes 7+(fe), J-{k) which determine the intensity of a<i and b_ 2 terms in a a (See 
Eq. fl3.3Q ). In Fig. |5.4j , we plot their dependence on momentum of emitted a particle. As we 



can see in Fig. |5.4| , 7+ amplitude falls down as k increases and goes to 2 as k — > 00. In this 
limit, 7+ defines mixing due to a simple rotation between a± and a 2 states. Since it can 
be successfully computed within the framework of quantum mechanics, it gives the usual 
Pontecorvo formula with only one oscillatory term. On the other hand, 7_ appears with an 
antiparticle creation operator and describe Bogoliubov rotation between a\ and 6_i states. 
This term is also responsible for ei ~^ 62 high frequency term and antiparticle beam creation. 
As we see in Fig. |5.4| , it decreases as k — > 00 and the mass difference becomes washed out by 
the relativistic gain of mass. This also means that at ultrarelativistic limit the QFT-mixing 
effects vanish so that the simple Pontecorvo formula is restored for flavor-oscillation. 
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FIG. 5.5. population density evolution for k=0.5GeV 
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To demonstrate the energy dependence, we show in Fig. p.5| the plot of population densi- 
ties evolving with time for the larger momentum k = 0.5 GeV. The line assignments are same 
as shown in Fig.|5.2|. As easily seen in Fig.|5.5|, the intensity of antiparticle beam decreases 



dramatically to about 10% (in contrast to 20-40% in Fig. |5.2j ) of initial intensity. The initial 
increase in the population density fluctuation in particle beams also reduces even though the 
quantum mechanical simple oscillations with ei ~ 62 frequency are still visibly distorted. Two 
beams nevertheless demonstrate strong correlation of the same kind as correlation in quan- 
tum mechanical EPR-paradox problem so that total flavor charge is preserved as it should 
be. It is also noticeable that there exist moments of time when the antiparticle outnumbers 



the particle of the same sort thus producing a negative flavor charge as shown in Fig. 5.3 



V. CONCLUSIONS AND DISCUSSIONS 



The non-trivial scalar and pseudoscalar meson mixing effects may be understood by 
the condensation of corresponding flavor states in the vacuum ||. We have now extended 
the analysis replacing the perturbative vacuum to the nonperturbative (flavor) vacuum. 
Central to this analysis is the interplay between the base (unmixed) Fock space and the 
physical Fock space. Their nontrivial relationship gives rise to the mixing and oscillation 
phenomena. While the similar quantum field theoretic formulation was presented for the 



fermion mixing P,lQ[], as well as boson mixing |]TT|, our analysis differs in the derivation of the 
normalization factor Z given by Eq.( [2.25| ) which is crucial to show the unitary inequivalence 
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between the mass vacuum and the flavor vacuum. We presented a new algebraic method 
which is distinct from the conventional method of using a differential equation for Z. While 
the unitary inequivalence occurs only in the infinite volume limit even for the boson case as 
discussed in Ref. ||11|| , we find an intrinsic difference between the fermion and boson cases. 
As shown in this work, the normalization factor Z for the boson given by Eq. fl2.25|) has 
a singularity on the complex plane at 7_ sin 9 = 2i while the corresponding result for the 
fermion doesn't have any singularity. As we summarized in Appendix B, this singularity 
corresponds to the divergence of the Taylor series expansion in powers of 9 for the region 
7_# > 2. For both the boson and fermion cases, however, the non-trivial observable mixing 
phenomena cannot occur unless there is both a nonzero mixing angle and also a nonzero mass 
(energy) difference between the two physically measurable mixed states. Dramatic small 
oscillations due to the virtual pair creation occur in the exact vacuum analysis, while only a 
simple sinusoidal Pontecorvo oscillations occur in the perturbative vacuum treatment. Some 
clarifying remarks on the Green function method and the arbitrary mass parametrization 



discussed in the previous literatures [|iy,|Il],|13|] are also summarized in Appendix C. 

As a physical application, we used our formulation to analyze the rj rf system and found 
that the measured mixing angle and mass difference between r] and 7/ can be related to 
the non-trivial flavor condensation in the vacuum. However, more fundamental questions 
such as the translation of the condensation in hadronic degrees of freedom to those in 
quark and gluon degrees of freedom remains unanswered. The answer to this question 
depends on the dynamics responsible for the confinement of quark and gluon degrees of 
freedom and perhaps has to rely on lattice QCD and/or some phenomeno logical model that 
accommodates strongly interacting QCD. Further investigation along this line is underway. 
Also, it would be interesting to look at the mixing transformations between gauge vector 
bosons governed by the Weinberg angle in the electroweak theory as well as vector mesons 
such as the p and uj. While the statistics are the same as the scalar and pseudoscalar bosons 
considered here, there will be additional spin dependent interactions which complicate the 
analysis. 
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APPENDIX A: COUPLED EQUATIONS FOR COMPUTING THE FLAVOR 

VACUUM STRUCTURE 



In this Appendix, we summarize the procedure for deriving Eq. fl2.17|) that describes a 



structure of the flavor vacuum. We first define flavor vacuum in terms of mass-eigenstates 
as the most general linear superposition of the form: 

\0;9>= Yl C nlmk (B)(a\) n (al)\bU) m (bU) k \0> 

n,l,m,k 

= E %f^(4r(4)'(& t - 1 ) m (^ 2 ) fe |o>. (A.i) 

n,l,m,k 

Then, using the definition of flavor vacuum: 

a a)/ g|0; 6 >= 0; 

b- a ,- p \0;9>=0, (A.2) 



and explicit expression for the ladder operators given by Eq.( |3.3| ), one can derive an infinite 
set of linear equations for C' n ik m coefficients: 



CC ?i,!+l,m)f — S + C n +l,lmk + S—C r d,rn-l,k = 0, 

n, l,m,k — 0, 1, 2, 3, . . . 

(A.3) 

To solve this infinite set of equations, we can express C'^ymt, C'n,J+i,fcm in terms of 
C'ni,m~i,k, C'nim,k-i so that we step by step reduce the n + I number of particles. Denoting 

(l:t:)- 5 -(-lc + r (a - 4) 

(notation for Z is chosen in correspondence to the index of particle type), we can write this 
relationship as 

C n+l,lmk = Z12C n lm,k-l + Z\\C nl,m-l,k 

C n ,l+l,mk = Z22C ' n lm,k-l + -^2lC"?ii,m-l,fc- (A-5) 

One also can write this in a symbolic manner introducing a kind of shifting operators with 
definition 

nlmk ^ nlm,k—lt 
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TTlC nlmk — C nl,m-l,k- (A-6) 

With the use of Eq.(A.6), Eq.(A.5) may be rewritten as 



C'n+i,imk — {Znk + Ziim)C' n i m k — (Zi 2 k + Z n rh) C"„_i 



Imk 



C'n,l+l,mk — (Z 22 k + Z 2 irh)C' n i m k — (Z 22 k + Z-xymfG 1 ' n ,l-l,mk — ■ ■ ■ > (A-7) 

and finally it can be written as 



n I 



C nimk — (Zi 2 k + Z u m) n (Z 2 2k + Z 21 m) l C' 00m k 



( E E Cn'Cf'Zu ^T m '4i^^" + ^ (m ' +t,) ^ m ' +t ')CW fe . (A.8) 

m'=0t'=0 

One should note that, since total momentum of vacuum state should be zero, C'oomk — 
unless m = k = 0. Therefore, in Eg. ( |A.8| ) only terms with {m! + t' = m,n + l — {m! +t') = k) 
must survive and from Eq. flA.4Tj we get Z\\ = —Z 22 , Z\ 2 = Z 2 i to find: 

oo n+l r5 

|0, 9 >= Z J2 E - J ^(S t i) n (4) Z (&li) m (6 + - 2 ) n+ ^ w iO >, (A.9) 



n,l=0m=0 n - 1 - 



where 



B nlm = E C^C t ;z[X m '- t 'Z^ ml+t \-l) l - t '. (A.10) 
m' + t' = m 
< m! < I 
< H < n 

Using a direct expansion, one can also verify that the above expression for vacuum state is 
equivalent to 

|0, 9 >= Zj2^(ZiAbU + Z 12 albl 2 ) n (-Z u albl 2 + Z 12 4bl 1 ) l \0 > . (A.ll) 



APPENDIX B: PERTURB ATI VE EXPANSION IN 9 FOR THE FLAVOR 

VACUUM 

In this Appendix, we try to directly estimate the norm of the flavor vacuum state 
G- 1 (9) |0 >= exp ( — 9S^j |0 > using the perturbative expansion in powers of 9 and show 
that the perturbative calculation of the flavor vacuum state is indeed impossible for large 
7_#. Truncating the series for G^ 1 (9) |0 > to the N terms, we have the term with the largest 
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number of particle coming from -y (a\b^_ 2 + a\ bl 2 ) in the (-9S) . Thus, the truncated 
series of G^ 1 (9) |0 > can be written as 



GJ, 1 (9) |0 >= X + i (-2=2 ) N (a\bl 2 + a\bU) N |0 > 



(B.l) 



= X + 7*T (-¥)" El» c» («■!)" (it,)" (4)""" lo >, 

where X denotes all terms with the total number of particles and antiparticles less than 2N. 
For the norm of above expression we can write then: 

'' " " ' " V . n\n] (N - n)\ (N - n)\ , v2 



G- N 1 (9)\0>\\ =||X|| 2 + (^) ^^„!„!(jy-n)!(N-n)!^f^ 

= \\X\\ 2 + (N + 1) (Y)'" 



Thus when 7_# > 2 the norm of the |O,0 >at= G^ 1 (0) |0 > is growing as the number of 
terms kept in expansion of the G (9) grows and therefore the transformation operator G^ 1 (9) 
is not well defined operator in the mass-eigenstates Fock space. 

One may also try to check directly the identity G (9) G^ 1 (9) — 1. In this type of approach 
one defines 



A' 



Mr 

G(9) = lim G N (9) = lim V ^ 

n=0 



Then, one shall prove that hin^^oo G^ (9) Gj/ (9) — 1 



0, 



(B-3) 



i.e. 



lim 



N^oo 



(G N (9) G^ 1 (9) 



\x > 



for any mass-eigenstate state \x > if G (9) is well 
defined. When multiplying Gn (9) and Gjf 1 (9) one typically get all coefficients vanished till 
the power of N and then have a "tail" up to the S 2N coefficient. If G (9) is well defined, 
this tail is expected to vanish when N is taken to infinity. However, this does not always 
happen in the perturbative expansion. To demonstrate this one may consider the last term 
of the "tail" given exactly by Jr^y. Recalling that S generator contains ^f- (a\b^_ 2 + a\b-i) 
combination, we can write the state (GnG^ 1 — lj |0 > as 

{G N G- K ' - 1) |0 >= Y + ^|^ («J)' (6»_ 2 )' (4) 2 "- (^f- |0 >, (B.4) 

where Y denotes all states with less then AN number of particles and antiparticles. The 
norm of this state is then given by 

2 ' V Vr >v (GL r ) 2 t\t\C2N n'r> V - /«,- 

(B.5) 



(GWGtf -l)|0>|| = \\Y\\ 2 + ^ (^) Y? t ^Ci N ft\t\(2N-t)\(2N-t)\ 

(2N + 1)(^) 2N . 



> 
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Again when > 2 the above expression is not convergent and the G (9) G^ 1 (9) expression 
is in fact not well defined in terms of mass-eigenstate fields. 

For small values of 9, however, the perturbative expansions are indeed convergent and 
the radius of convergence is related to the pole of Z =< 0\0,6 > on the complex plane of 9. 
Here, the pole (critical) value is given by 7_ smh(9 cr iu C ai) = 2. 

APPENDIX C: REMARKS ON THE GREEN FUNCTION METHOD AND THE 
ARBITRARY MASS PARAMETRIZATION 

1. Green Function Method 

We note that a straightforward use of the Green function with the conventional definition 
< O|T[^(x) , 0(y)]|O > encounters some difficulties in the mixing analysis due to the fact that 
the flavor vacuum state is not stationary in time (|0, 9 > (t) ^ |0, 9 > (£'))• The conventional 
Green function cannot be adopted without specifying at which times the flavor vacua were 
taken in the inner product. In fact, the most obvious generalization of the Green function 
as the overlap between the states created at times x° and y° (i.e. G(a — > /3; x°, y°) 
=< 0, y°\T[^/3xjj a ] |0, x° >) breaks down due to the unitary inequivalence of the flavor Fock 
spaces at different times. Therefore, some sort of modification, like parallel translation 
of states to the same time, shall be needed to define the Green function appropriate in 
the mixing analysis. The flavor mixing problem can then be treated using this modified 
propagation functions as discussed in the previous literature [ |T0|| . 

In the process of our calculations we also noticed that some entities indeed appeared as 
"transition" amplitudes from one state to another. For example, < N a >= Z a + |A*a| 2 
can be considered as superposition of "vacuum rotation" background contribution Z a and 
contribution from a — > a transition with A aa transition amplitude. In this manner, one can 
introduce the Green function that only accounts for the transition amplitude without the 
vacuum contribution. In this way Green function can be defined by: 

G aa (x, t; y, 0) =< 0>„(z, t)^ a (y, 0)|0' > (C.l) 

where the vacuum state is taken at any (but certain) fixed time, for example t = 0, and this 
coincides with the definition given in ||10|| . For the propagator with a definite momentum k, 
we then obtain: 

G aa (k,t)=A* aa (k,t), 
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G aP (k,t)=A* pa (k,t),etc. (C.2) 

Eq. ( |3.7| ) then allows to define the propagation functions for any kind of transition. 

We should note however that the treatment with such a modified Green function does 
not cover all the variety of the effects in the mixing problem. In particular the condensate 
contribution Z a , that can be related to the unitary inequivalence of the flavor Fock spaces 
at different times, is lost so that this part of problem is missing when the above approach is 
taken. Nevertheless, the Green function method is useful in calculations of flavor-operator 
expectation values, scattering amplitudes and other quantities in which vacuum contribution 
Z cancels out. 



2. Arbitrary Mass Parametrization 



In this sub-Appendix, we remark on the arbitrary mass parametrization flll"l , |T3"ll i n the 



mixing problem. As discussed in Ref. ||13|| , one may treat the flavor fields that were initially 
written as 

dk 



(2tt) 3 / 2 

equally well in an arbitrary mass basis, i. e 

dk 



(C.3) 



(u kji a k:i (t) + £_ feji 6i fcji (t)) 



Jkx 



(C.4) 



(2tt) 3 / 2 

where u k ,i and v k ^ are free-field amplitudes with some new arbitrary masses. Since there is 
no physical reason to prefer one form over the other, Ref. [T3| and then Ref. [IT] claimed 
that no arbitrary mass parameters should appear in physically observable quantities, i.e. 
they shall be invariant under specific Bogoliubov transformation going from Eq.( |C.3| ) to 
Eq.( 
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(C.5) 



It is true [|Tl~1 , |l3f that the perturbative vacuum treatment |||| yields the normalization 
of the flavor state not as unity but as some constant that depends on the arbitrary mass 
parameter. In particular, Eq. (|3.3j ) can be viewed as an expansion of the flavor ladder operator 
in some basis constructed from the free-field ladder operators. Then, the normalization of the 
one-particle state in perturbative vacuum treatment |9j was given only by IS*-! 2 coefficient 
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at the b[ 2 operator, that is obviously depending on the choice of basis, e.g. changes with 
rotation in {(ai, a2)&(6i, 62)} plane. Such arbitrariness is completely avoided in the exact 
vacuum treatment because the normalization of the flavor state is given by unity no matter 
what mass-basis is used. 

However, the claim [II],[n| that the number expectation values are not physical because 
they do depend on the arbitrary mass parameters cannot be correct. As a counter example 
to such claim, one can consider a very specific case of the mixing problem, namely when the 
mixing is absent (G(9) = 1). As we discuss below, applying such claim ]TT| , |I3| to this specific 
example leads to a conclusion that cannot be correct. With no mixing, we are dealing with 
nothing else but free- field problem where the particle number operator Ni(t) = a\{t)<ii(t) and 
the particle number expectation value < iVj >= |{aj(0), a|(t)}| 2 (for simplicity we consider 
fermion case) must be well-defined physically observable quantities. However, when we 
apply the transformation given by Eq. (|C.5|) and compare directly Eq. (|C.3|) and Eq.( |C.4j ), 
we observe that 

5< = (4,^MKie~" M + (^-U-k^bU/^ = Pka^e-^ 1 + Xkbl,/ 6 ^ (C.6) 

and 

< >= |{5,,aJ(t)}| 2 = IWl'e-^ + |A*|V^f , (C.7) 



where we follow the notations in Refs. |Tl|J13[| . In the case of free fields the number expec- 
tation value does depend on the arbitrary mass parameters and thus following the above 
claim one may conclude that < iVj > is not a physically measurable quantity. However, 
this cannot be correct because both the particle number operator and the number expec- 
tation value in the free-field problem are well-defined physical observables. We viewed the 
above inconsistency as follows. The transformation given by Eq.( |C.5| ) is in fact nothing else 
but redefinition of the particle states, so that the tilde quantities correspond to some new 
quasiparticle objects and the number operator now describes the number of different type of 
particles than before. Therefore, the number operator average shall not be expected to be 
same in such transformations. Indeed, it should change in some covariant and self-consistent 
manner instead of being invariant under such a redefinition. 
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